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ABSTRACT 

The number count of ionized bubbles on a map of 21 cm fluctuations with 
the primordial non-Gaussianity is investigated. The existence of the primordial non- 
Gaussianity modifies the reionization process, because the formation of collapsed ob- 
jects, which could be the source of reionization photons, is affected by the primordial 
non-Gaussianity. In this paper, the abundance of ionized bubbles is calculated by us- 
ing a simple analytic model with the local type of the primordial non-Gaussianity, 
which is parameterized by /jvl- In order to take into account the dependence of the 
number count on the size of ionized bubbles and the resolution of the observation 
instrument, a threshold parameter B\, which is related to the the surface brightness 
temperature contrast of an ionized bubble is introduced. We show the potential to put 
the constraint on /jvx from the number count by future observations such as LOFAR 
and SKA. 



1 INTRODUCTION 



Recent cosmic microwave background (CMB) observations 
have revealed the statistical nature of primordial fluctua- 
tions a nd given the strong s upport of the inflationary sce- 
nario (|Komatsu et al.l I2009T ) . However, it is still difficult 
to specify the model among the inflationary scenario from 
CMB observations. The measurement (or even the upper 
bound) of non-Gaussianity of primordial fluctuations is, on 
the other hand, expected to have the potential to rule out 
many inflationary models. Although the primordial fluctu- 
ations are predicted as nearly Gaussian in all inflationary 
models, the degree of the deviation from the Gaussianity de- 
pends on each specific model. The slow-roll inflation models 
with a si ngle scalar field p r oduce almost pure Gaussian fluc- 
tuations (|Guth fc Pilll982l ; IStarobinskvlll982l ; | Bard een et al.l 
1983), and the deviation from the Gaussianity is unobserv- 
ably small jFalk et al.lll993l ; lGangui et alJll994 ). Meanwhile 
some multi-field inflation scenarios can produce large non- 
Gaussianity which can be observed in the near future ex- 
perim ents such as PLANCK mission (fBattcfcl d fe Eastherl 
120071 ). For comprehensive review see lBartolo et alT "( 2004 ). 



The most stringent constraint on the primordial non- 
Gaussianity is cur rently obtained from CMB measurements 
(| Verde et al.ll2000h . Wilkinson Microwave Anisotropy Probe 
(WMAP) sets the limit on the so-called local type of the 
primordial non-Gaussianity, which is parameterized by the 
constant dimensionless parameter /jvl- From WMAP 5yr 



data, the constraint was —9 < /jvl < 111 (Komatsu et al 
2009). Recently, this constraint is updated by ICurto et al 
(2009), -18 < /jvl < 80. 



The large scale structure and the abundance of the col- 
lapsed objects, e.g. galaxy clusters and galaxies, are alter- 
native probes of the non-Gaussianity. The positive /jvl en- 
hances the abundance of the collapse objects, while the neg- 
ative /jvl decreases the abundance. In particular, the for- 
mation process of rare objects such as high-mass collapsed 
objects at high redshift, which is contr olled by the high den- 
sity t ail, is quite sensitive to /jvl (see iDesiacques fc Seliakl 
l20ld for a review). On the contrary, the abundance of a 
void, which is formed in the underdense region, is sensitive 
to the low density tail of the fluctuation distribution, and has 
also the potential t o prove the primordial no n-Gaussianity 
(|Grossi et al.ll2008l ; iKamionkowski et al.ll2009l ). 



The non-Gaussianity also affects the reionization pro- 
cess, because the formation of high-mass collapsed objects, 
which could be the source of reionization pho tons, is en- 
hanc e d (reduced) by posit ive (negative) /jvl (|Chen et alj 
120031 ; I Crociani et al.1 |2009| ). Therefore, in this paper, we 
study the potential of the number count of ionized bub- 
bles on a map of the redshifted 21 cm line fluctuations 
from neutral hydrogen as a probe of the non-Gaussianity. 
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Now, LOFAK0, and SKA0 are being installed or de- 

signed for the measurements of 21 cm line fluctuations. The 
maps of 21 cm fluctuations are sensitive to the density, tem- 
perature, and ionized fraction of the intergalactic medium 
(IGM). Studying the 21-cm tomography tells us about the 



epoch of reionization. iMadau et al.l 1997 


; Tozzi et al.l 


2000; 


ICiardi & Madaul 120031; Furlanetto et all 


20041). An ionized 



bubble will be observed as a dark spot on a map of 21 cm 
fluctuations, because hydrogens are fully ionized in the in- 
side of the bubble. We employ a simple analytic model of 
the ionized bubble to investigate the non-Gaussian effect on 
the number count of ionized bubbles. 

The outline of this paper is the following. In Sec. II, we 
give the simple analytic model of the ionized bubble based on 
the Press-Schechter theory with the non-Gaussianity /jvl- 
In Sec. Ill, we calculate the surface brightness temperature 
of a ionized bubble and show the results of the number count 
of ionized bubbles. Section IV is devoted to the conclusion. 
Throughout the paper, we use the concordance cosmological 
parameters for a flat cosmological model, i.e. h = 0.73 (Ho = 
h x 100 km/s/Mpc), T = 2.725K, fi b = 0.05, fi m = 0.27 
and erg = 0.8. 



2 ANALYTIC MODEL OF IONIZED BUBBLES 

In this paper, we focus on the so-called local type of the pri- 
mordial non-Gaussianity, which is para meterized by the con- 
stant dimensionless parameter /jvl as jKomatsu fc Spergell 
l200ll ) 



$(X) = $ G (X) + f!fL(Q%(x) - {$%(x))), 



(1) 



where $ is Bardeeen's gauge-invariant potential, $g is the 
Gaussian part of the potential and ( ) denotes the ensemble 
average. The power spectrum of <&g is defined by 



(<E> G (fc)$ G (fc')} = (2ir) 3 S D (k + k')P(k), 



(2) 



where 8d is Dirac's delta function. The bispectrum of the 
local type potential can be written as 

($(fei)$(fe 2 )*(fc 3 )> = (27r) 3 5 D (fe 1 + k 2 + fe 3 )B(fei, fe 2 , fcs), 

(3) 

where B(ki, k 2 ,k 3 ) can be written with /jvl and the power 
spectrum P(k), 

S(fei,fe 2 ,fe 3 ) = 2f NL (P(k 1 )P(k 2 )+P(k 1 )P(k 3 )+P(k 2 )P(k 3 )). 

(4) 

The existence of the primordial non-Gaussianity modi- 
fies the mass function obtained under the assumption of the 
Gaussian density fluctuation. The mass function for the non- 
Gaussian case is given by estimating the fractional correc- 
tion to the one for the Gaussian case in the Press-S chechter 
theory (|Matarrese et al.ll2000l ; I"Lo Verde et alj|2008h . 
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2 http://www.mwatclescope.org/ 
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where IZng is the non-Gaussian c orrection. We adopt the 
correction of lLo Verde et all l|2008l ). 
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where S c is the critical density contrast for collapsing ob- 
jects, o m and S3 are the smoothed density dispersion and 
the smoothed skewness at a redshift 2 with a top-hat window 
function W(R, k). The smoothed density dispersion <jm(M) 
is given by 



a 2 M (M) = (4) 



dk A 



■W(R,kyD(k,z) 2 P(k), (7) 



where R is the scale corresponding to M in the top-hat win- 
dow function and D(k, z) is the relation between $ and 
8 in the linear perturbation, 8(k,z) = D(k,z)&(k). The 
smoothed skewness Ss(M) is obtained by 



Ss(M) 

<4> 
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(S 2 R ) 2 ' 
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(8) 



(2.)3(2.)3(2.)3 W/ ^ fcl ^^ fc2 ^^ fc3 ) 
xD(k 1 ,z)D(k 2 , z)D(k 3 , «)(*(fci)*(fc)*(*8)).(9) 



The modification of the mass function from the existence of 
non-Gaussian fluctuations is described by the terms propor- 
tional to S3. 

The reionization involves various physical processes 
such as star formation, radiative transfer through IGM, and 
so on. The aim of this paper is to demonstrate the poten- 
tial of the number count of ionized bubbles as a probe of 
the non-Gaussianity. Without discussing the details of the 
reionization process, we employ a simple analytic reioniza- 
tion model. First, we assume that a halo whose virial tem- 
perature is larger than 10 4 K can collapse and the source 
of ionizing photons can be produced in the inside of such 
collapsed halos. Secondly, the shape of the ionized bubble 
produced by a collapsed halo is spherical and the size de- 
pends on the halo mass. Under these assumptions, we give 
the ionizing photons per baryon in a collapsed halo as 



iV 1( 



(10) 



where N 1 is the number of ionizing photons per baryon in 
stars, /„ is the star formation efficiency, and / eS c is the es- 
cape fraction of ionizing photons into the IGM. We adopt 
/» = 0.1 and / osc = 0.05 in this paper. The parameter JV 7 
depends on the ionizing source. Here we consider popula- 
tion III stars as the ionizing sources and take JV 7 = 44, 000 
ijBromm et al.l l200lT l. Assuming that all ionizing photons 
ionize hydrogen atoms and recombinations are negligible, 
we obtain the maximum physica l size of an ionized bubble 
in the IGM as ijLoeb et al.ll2005T ) 
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Figure 1. The sketches of ionized bubbles and the smearing ef- 
fect by the finite observation resolution. Ionized bubbles are fully 
ionized in the inside of each radius R(M) (top panel). The finite 
resolution smoothes the signal from ionized bubble. The smoothed 
surface brightness temperature contrast B for each bubble is dif- 
ferent, depending on both the size of ionized bubbles and the 
resolution size (bottom panel). As a result, each threshold Bb 
gives different number count of ionized bubbles. 



3 NUMBER COUNT OF IONIZED BUBBLES 

An ionized bubble is detected as a hole with radius f? max 
on a map of 21 cm fluctuations because the inside of an 
ionized bubble is assumed to be fully ionized. In an actual 
observation, however the hole is smeared due to the finite 
beam size of a telescope with the signal from surrounding 
neutral hydrogens. Taking into account this smearing effect, 
we define the surface brightness temperature contrast of an 
ionized bubble as 

Y obs (6,z) 



B(0,z) 



T 21 {z)jdn> exp 



(12) 



where a = Ofwhm /V81n2 with #fwhm being the FWHM 
resolution of the observation, and T2i(z) is the 21 cm back- 
ground temperature at z. The surface brightness tempera- 
ture of an ionized bubble Y b s (9) is obtained by 



Y obs {6,z) = j dQ' T{9-6',z)exp 



(9-9'f 



2a 2 



(13) 



where the angular profile of the ionized region T(9, z) is 
given by 



T{e,z) = 



0, 9 < 9 R 
T 2 i{z), 9>9 R 



(14) 



Here 9r is the angular size corresponding to R ms , x - 

As the criterion for the detection of an ionized bubble, 
we introduce a parameter Bb. We count bubbles with B 
smaller than Bb as the detectable ones on a 21 cm map (see 
also Fig. [[J. Considering Eqs. (fTT|) . (Tj"3j) and (TH}, we can 
relate Bb to the limiting mass Mn m which is the mass of a 
collapsed halo associated to an ionized bubble with B = Bb- 
Therefore, the number count of ionized bubbles with B < Bt 
is given by 



N(< B b 



M lim (S b ) 



, , , dV dn . 
dM—-—Az, 

dz dM 



(15) 



where V is a comoving volume and Az is a width of the 
observation redshift slice. We set Az = 0.1 throughout this 
paper. 

In Fig. [21 we plot N(< Bb) as a function of Bb for 
different /jvl- The 21 cm fluctuations during the epoch of 
reionization will be mapped by radio interferometric obser- 
vations. The resolution #fwhm for a wave length A depends 
on its baseline D as #fwhm = X/D. In the calculation of 
Fig. [2] we choose two types of the baseline; low resolution 
type, D = 2 km, which corresponds to the LOFAR baseline, 
and high resolution type, D = 5 km, which corresponds to 
the SKA baseline. 

As is shown in Fig. [2] the number count strongly de- 
pends on the resolution and the redshift of the observa- 
tion. The large angular resolution smears the ionized bubble 
whose size is smaller than the resolution. Then, the more the 
redshift decreases, the larger the number of halos which grow 
enough to collapse is. In addition, both physical and angular 
size of a bubble at high redshifts is small compared to the 
one with the same halo mass at low redshifts. As a result, 
the number count strongly decreases with increasing of the 
angular resolution or the redshift. The number count with 
Bb — 0.1 for the high resolution type at z = 11 is about 
2000, while the number count for the low resolution type at 
z — 13 is less than one even for full sky survey with Bb = 0.5 
as is shown in Fig. [3] 

To make clear the difference of the number counts be- 
tween the Gaussian and non-Gaussian primordial fluctua- 
tions, we plot the ratio of the number counts between non- 
Gaussian and Gaussian cases as a function of Bb in Fig. [31 
Producing ionized bubbles during the epoch of reionization 
is a rare event. Therefore the number counts depend on the 
degree of the primordial non-Gaussianity. This tendency is 
emphasized as the redshift increases. The deviation from the 
Gaussian case becomes larger for all non- Gaussian cases at 
z — 13 than at z = 11. The difference is also enhanced 
with decreasing Bb because lower Bb corresponds to higher 
threshold halo mass Mu m . 

It is shown in Fig. [3] that the dependence on /jvl of the 
number count is more significant for the negative value than 
the positive one. In the case of negative /jvl, the high density 
tail of the distribution function is suppressed, and less num- 
ber of bubble formation, which is a rare event, takes place. 
Therefore with decreasing Bb, the number count rapidly de- 
creases as is shown for /jvl = —50. On the other hand, 
for positive /jvl, the high density tail of the distribution 
function is enhanced. Accordingly, there are more bubble 
formation for larger /jvl- However, the terms proportional 
to 5*3 and dSz/dM in Eq. ^ have opposite sign in the case 
of positive /jvl- Therefore the dependence on /jvl is not as 
strong as in the case of negative /jvl- 

By comparing upper two panels with lower two panels 
of Fig. [3] we find that deviation from the Gaussian case is 
larger in the lower resolution observation. However, number 
counts become smaller for the lower resolution case. As a 
result, the error of the determination on /jvl is expected to 
be large in the low resolution observation. We discuss this 
point in the term of signal to noise ratio of the detection 
later. 

In Fig. [J] we show the evolution of the deviation from 
the Gaussian case with fixing Bb- In these figures, we adopt 
the high resolution type and set Bb = 0.5 and 0.1. The 
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Figure 2. The number count of ionized bubbles as a function 
of Bf, for different /jvl- We plot the number counts for z = 11 
and z = 13 as the thick and thin lines, respectively. We use the 
low resolution type (D = 2 km) in the top panel and the high 
resolution type (D = 5 km) in the bottom. 

figures tell us that increasing redshift or decreasing Bb makes 
the deviation large. In particular, in high redshift or low Bb 
cases, the number count for the negative /jvl case damps 
rapidly due to the suppression of the bubble formation. 

In order to discuss the measurement of /jvl by the num- 
ber count of bubbles, we calculate the Signal to Noise ratio. 
We assume that the observation instrument is ideal, i.e., 
there is no instrumental noise, and the foreground noise can 
be completely removed. Therefore we take into account only 
the Poisson noise of the number count. The SN ratio for fixed 
Bb is given by 

S/N = ^fs ky N{< B b ), (16) 

where f 3 k y is the fraction of sky observed. According to 
Eq. ()16p . the large number count gives high SN ratio. There- 
fore, the larger the threshold Bb is set, the higher the SN 
ratio becomes. 

However, setting high threshold makes two problems 
arise. The first problem is difficulty counting bubbles. In the 
large Bb case, we need to count not only large bubbles but 
also small bubbles which are hardly seen in a 21 cm map. 
The second problem is that the deviation from the Gaus- 
sian case becomes smaller as Bb increases as is shown in 



Fig. [3] It should be notice that larger S/N does not nec- 
essarily mean sensitive to the non-Gaussianity detection. 
Therefore we need to optimize for the threshold value Bb 
corresponding to the observational settings. 

In Fig. we plot the SN ratio with adopting SKA's 
field of view (~ 100 deg 2 for 140 MHz), which corresponds 
to f sky = 2.42 x 10" 3 . The top panel of Fig. O is for the 
low resolution (D = 2 km) and the bottom one is for the 
high resolution (D — 5 km). In order to obtain S/N > 5 
for /jvl = 100 for the low resolution at z = 11 we must 
set B b > 0.7 at least. With B b ~ 0.7, according to Fig. [2l 
the number counts for /jvl = 100 is ~ 20 for SKA's field 
of view, while the one for /jvl = are ~ 10 for SKA's 
field of view (it should be reminded that Fig. [2] shows the 
number counts for the full sky survey, f a k y = 1). On the 
other hand, in the high resolution, Bb > 0.1 is enough to 
obtain S/N ~ 5 for /jvl = 100. In this case, the number 
count for /jvl = 100 is ~ 20, while the one for /jvl = is 
~ 10 with S/N ~ 3. Therefore we conclude that the number 
count of ionized bubbles has potential to give the constraint 
on /jvl • 



4 CONCLUSION 

In this paper, we have studied the number count of ionized 
bubbles on a 21 cm map as a probe of the primordial non- 
Gaussianity. Although we have used a simple analytic model 
of an ionized bubble, such a simple model is enough to il- 
lustrate the relation between /jvl and the number count of 
ionized bubbles. 

In order to take into account the dependence of the 
number count on the size of ionized bubbles and the res- 
olution of the observation instrument, we have introduced 
a threshold parameter Bb- The bubbles with B < Bb are 
counted as the detectable ionized bubbles on the 21 cm 
map. Setting smaller Bb means counting larger size bubbles. 
Therefore, we can expect large numbers of ionized bubbles 
for large Bb- The deviation of the number count for the non- 
Gaussian case from the one for the Gaussian case depends 
on the threshold Bb- Small Bb brings large deviation. How- 
ever the number count for small Bb becomes small. As a 
result, the SN ratio goes to low. In order to constrain /jvl 
by the number count of ionized bubbles, the optimization 
between Bb and the resolution of observations is needed. 
For example, in our model, we have found that the number 
counts for /jvl = 100 and /jvl = for the high resolution 
with Bb = 0.1 are about 20 with S/N ~ 5 and about 10 
with S/N ~ 3, respectively. 

The number count also depends on the redshift of a 
21 cm map. Compared with the number counts at low red- 
shifts, the ones at high redshifts are lower sensitive to /jvl. 
However the number counts at high redshifts are relatively 
smaller than at low redshifts. Therefore, the SN ratio at high 
redshifts becomes large and it becomes difficult to give the 
constraint on /jvl- Besides, the detection of ionized bubble 
itself is a challenge (|Datta et al.ll2007l , |2009| ). 

In this paper, we have not taken into account the effect 
of overlapping of ionized bubbles. The overlaps take place 
as the universe evolves. As is shown in many reionization 
simulation works, once overlapping proceeds, the number 
count of ionized bubbles quickly decreases and the size of 
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Figure 4. The evolution of the fractional deviation from the 
number count for the case of the Gaussian mass function as a 
function of redshift. We set Bj, = 0.5 in the top panel and B b = 
0.1 in the bottom panel. The curves are for /jvx = 100, /jvz, = 10, 
/jvl = —10 and /jvl = —50 from top to bottom in both panels. 



ionized bubbles becomes rapidly increasing. Moreover the 
evolution of the overlapping process depends strongly on the 
model of reionization. Therefore, our results in this paper 
can be only adapted for the early phase of the epoch of 
reionization before the overlapping process takes place. 

In order to properly treat the IGM ionization process 
including overlapping of bubbles, a detailed numerical sim- 
ulation with radiative transfer is needed. In the near future, 
we intend to perform the numerical simulation, considering 
the reionization process. 



0.01 



0.1 
B b 



1.0 



Figure 3. The fractional deviation from the number count for 
the case of the Gaussian mass function as a function of B b . We use 
the low resolution type (D = 2 km) and the high resolution type 
(D = 5 km) in the top two panels and the bottom two panels, 
respectively. In both angular resolution cases, we set z = 13 and 
2 = 11. The curves are for /jvl = 100, /jvl = 10, Jml = —10 
and /jvi = —50 from top to bottom in both panels. 
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